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ABSTRACT 


Steady-state harmonic stress waves in an isotropic elastic plate excited 
on one face by a circular transducer are analyzed theoretically. It is assumed 
that the transmitting transducer transforms an electrical voltage into a 


uniform normal stress 

at 

the 

top face 

of the plate. To 

solve the boundary 

value problem, first the 

radiation 

into a 

half- space is 

considered . 

This 

introduces longitudinal 

(P) 

and 

shear 

(S) waves into 

the plate . 

Then, 


reflections are considered successively at the bottom and top faces of the 

i 

plate . Each reflection produces both P and S waves for each incident P or S 
wave . 

A separate circular receiving transducer, also located at the top face 
of the plate, is considered. It is assumed that the receiving transducer 
produces an electrical voltage proportional to the average spatially 
integrated normal stress over its face due to an incident wave. 

The asymptotic behavior of the frequency response at a receiving point 
due to a multiply reflected wave is given. It is found that a receiving point 
observes an incident wave asymptotically as a plane wave which propogates and 
reflects in the direction of the multiply reflected ray constructed 
geometrically using Snell's law. The far field condition, for which the 
asymptotic solution is valid, is also discussed. This condition suggests that 
although the thickness of the plate may be small, the waves which reflect 
sufficiently many times do satisfy the far field condition. 

A numerical procedure is given to evaluate the frequency response at a 
receiving point due to a multiply reflected wave in the near field. Its 
stability and convergence are discussed. Also, exponential decay is introduced 



IV 


to account for material attenuation. 

Calculations are done for aluminum plates. It is found that the 
numerical procedure becomes unstable for h/X^<1.6 where h is the thickness of 
the plate and X^ is the P wavelength. 

Parametrized plots which determine the particular wave whose frequency 
response has maximum magnitude compared with other multiply reflected waves 
are given for a range of values of dimensionless parameters. The effects of 
changes in the values of the parameters are discussed. 
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1 INTRODUCTION 


An ultrasonic testing (UT) parameter called the stress wave factor (SWF) 
was recently introduced fay Vary et al. Cl , 21 for nondestructive evaluation 
(NDE) of materials. In the SWF configuration, separate transmitting and 
receiving transducers are coupled to the same face of the test structure and 
the number of oscillations exceeding a preselected voltage threshold (or some 
modifications thereof C33) in the output signal, due to an input pulse having 
a broadband frequency spectrum, is defined as the stress wave factor. Williams 
and Lampert C33 indicated a correlation between the SWF and the attenuation in 
carbon fiber composites by showing that as the residual strength decreases due 
to impact damage, the SWF decreases and the through- transmission ultrasonic 
attenuation increases. Williams et al. have also shown that the initial 
through- transmission longitudinal attenuation can be correlated with the 
compression fatigue life C4] and the flexural fatigue life [53 of carbon fiber 
composites . 

Unlike conventional pulse- echo testing C63 where nonoverlapping 
reflected wave echoes are analyzed, the SWF is also valid for the analysis of 
overlapping echoes. 

Although the SWF has been generally used to characterize microstructural 
defect states of materials, it may also be used to detect cracks and 
delaminations fay coupling separate transmitting and receiving transducers to 
the same face of the test structure. The output in such a SWF configuration 
may be significantly affected by cracks oriented perpendicular to the faces of 
the structure, which are otherwise difficult to detect by most NDE methods. 

An important step in any quantitative ultrasonic NDE procedure is the 
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analysis of the stress wave transmission characteristics of the test structure 
through which the transmitting and receiving transducers communicate. The 
ultrasonic input- output characteristics of the SWF configuration containing 
thick isotropic elastic plates are studied in Z71 theoretically and 
experimentally . It is assumed that the transmitting transducer transforms an 
electrical voltage into a uniform normal stress on the plate and vice versa 
for the receiving transducer. The asymptotic value of the normal stress is 
calculated for an isotropic elastic half- space subjected to a uniform 
harmonic normal stress applied to a circular region at its surface, using the 
results of Miller and Fursey C83 . The top surface of the plate is assumed to 
coincide with the surface of the half— space and the bottom surface of the 
plate is assumed to be at a depth within the half— space equal to the plate 
thickness. Then the radiated stress waves are traced within the plate by 
considering reflections at its top and bottom faces. Reflection coefficients 
for plane waves are used assuming that the thickness of the plate is so large 
that the bottom face of the plate is in the far field of the transmitting 
transducer and the asymptotic spherical waves are nearly plane at the point of 
reflection . 

The purpose of this study is to extend the analysis for thick plates 
given in C7] to plates of any thicknesses; in particular, to plates whose 
thicknesses are smaller than the far field distance of the transmitting 
transducer . 

One way to address the boundary conditions which are imposed by the 
transmitting transducer at the top face and the boundary conditions at the 
stress— free bottom face is to expand the radially propagating plate waves 
which are derived by the nonlinear Rayleigh- Lamb frequency spectrum C9] . Such 
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an approach would be especially useful at low frequencies where the thickness 
of the plate is of the order of the longitudinal wavelength because a small 
number of radially propagating waves would be excited. 

Another way to approach the problem, which is more practical at high 
frequencies where the thickness of the plate is much larger than the 
longitudinal wavelength, is to consider the radiation from the transmitting 
transducer into a half- space and then to introduce succesive multiply 
reflected waves at the bottom and top faces of the plate. Each multiply 
reflected wave whose potential is expressed in integral form is called a 
generalized ray. This method is often called the generalized ray theory and is 
discussed in C103. Transient waves generated by a variety of internal and 
surface forces in a plate have been analyzed by Pao and Gajewski C103, and 
Ceranoglu and Pao Cl 11 using this approach. 
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2 STRESS WAVES IN THE PLATE 


2 . 1 Governing Equations and Boundary Conditions 

Consider an isotropic elastic plate of thickness h, bounded by the z = 0 
and z=h planes and extending infinitely in the x and y directions as shown in 
Fig . 1 <a) . As done in C73 , it is assumed that the transmitting transducer 
transforms an electrical voltage into a uniform normal stress on the plate. 
The cartesian coordinates x,y,z; the cylindrical coordinates r,#,z; and the 
spherical coordinates R,©,# are shown in Fig. 1 (b) . 

Then, the boundary conditions on the top face of the plate, including 
those imposed by the transmitting transducer, for steady-state conditions can 
be written as 


at z = 0 


o' 

z z 


e 

for 

0< r< a 

0 

for 

C>a i 


< 1 > 


o' - or 
r z 


0 


where o' , o' / O' are the components of the stress tensor in the cylindrical 

coordinates C123 , w is radian frequency, t is time, -1 , and a^ is the 

radius of the transmitting transducer. The amplitude of o' is taken as the 

zz 

unity since only the frequency response of the plate is of interest. Because 
the bottom face of the plate is stress- free, the boundary conditions there 


can be stated as 
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at z=h : 


o' = o' = cr 
zz rz z# 


0 


< 2 ) 


Because of the circularity of the excitation/ there is axial symmetry 
about the z axis and so the displacement component in the & direction 
vanishes. Then, the theory of wave propogation suggests that the following 
uncoupled wave equations are satisfied within an isotropic elastic solid 1121 


1 3 2 ® 



( 3 ) 


V 2 <p 


1 3 2 *P 



(4 ) 


where ® is the Iaplacian operator and for cylindrical coordinates with axial 
symmetry is given by 


= 


1 3 3 4 


3 r 2 r 3 r 3z 2 


( 5 ) 


<t> and ® are uncoupled longitudinal <P> and shear (S) wave potentials; and c 
and c 2 are the P and S wave velocities, respectively. The components of the 
displacement and stress fields can be found from the potentials and the 
relations are summarized in Table 1 . 

The complete solution requires finding the wave potentials which satisfy 
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either eqn. (3) or eqn. (4) such that the superpositions of the stresses 
generated from these potentials via the relations given in Table 1 satisfy the 
boundary conditions at z=0 and z=h as given by eqns . <1> and (2), 

respectively . 


2.2 Half-space Solution 


At the first stage of the solution, the boundary conditions at z=h are 
ignored. Then the problem reduces to that of radiation into a half- space. 
Using Hankel transform techniques, Miller and Pursey C8 ] solved this problem 
for the boundary conditions given in eqn. <1) . That solution is reviewed in 
C131 with the notation used in this study and is given by 


P 


roo r z . 

= — I E b (I ) J n (2TT5 ) exp ( 2fi i J 1 ~S 2 ) dS 

1 J Q 


( 6 ) 


X, 


and 


2tt a /*w 


c b = 


S 


l± f 

1 •'O 


r X r 

,<S> J ( 2*tt 5 ) esp(2TTi 


2 2 

k -5 > d€ 


< 7 ) 


where <t> p and are the P and S wave potentials which satisfy eqns. <3) and 
(4), respectively. The shear modulus of the plate material is X is the P 
wavelength and is given by X^ =2itc j/«; € is a dimensionless integration 

variable introduced in the Hankel transform; J Q denotes a Bessel function of 
the first kind and of order zero; k^c^c^ E p (S) and ) are called the 

excitation functions for the P and S waves, respectively, and are given by 
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E p (€ > 


2 2 

C2€ -k ) 


G < € > 


J ( 2 ITS ) 

X 1 


E s (€ } = “ 2 1 


J 1-S 2 a 

J (2ttS ) 

G <S ) 


( 8 ) 


( 9 ) 


where 


G(5 ) = 


<2S 2 -k 2 ) 2 


4€ 


2 J l-S 2 J k 2 -s 2 


( 10 ) 


and J is a Bessel function of the first kind and of order one. 

The potentials given in eqns. (6) and (7) and the potentials which are 
introduced below all have a steady- state harmonic time dependence, indicated 
by exp(-iut); this harmonic time dependence function is dropped from the 


expressions for convenience . 

The excitation functions in eqns. <8) and (9) have a pole at € =€ where 

the denominator function, G<€ ) , given in eqn. <10) vanishes when 5 R * c R / c 1 

where c is the Rayleigh wave velocity. € is the so- called Rayleigh pole. 

R ** 


2.3 Reflected Waves 

The superpositions of the stresses generated from the potentials <t> p and 
satisfy the boundary conditions at z=0 only. So, in the so- called first 
stage the boundary conditions at z=h are violated. In the second stage, new 

waves which are called reflected waves are introduced so that the boundary 

conditions at z=h are then satisfied. This requires introducing a set of P and 
S waves for each P and each S wave initiated in the first stage, which is due 

to mode conversion [141. These new waves are called the PP and PS, and SP and 

SS waves and their potentials are denoted by <t> pp and <t> ps , and <t> sp and <t> ss , 
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respectively. In particular the potentials ^pp and o are introduced to 

cancel the stresses o' and cr at z=h generated from the potential O , 

zz rz p 

corresponding to the reflection of the P wave at the stress- free plane 
boundary at z=h. 

The superpositions of the stresses generated from the potentials 
introduced in the first and second stages satisfy the boundary conditions at 
z=h only. In the third stage, additional reflected waves are introduced so 
that the boundary conditions at z=0 are then satisfied. Again, this requires 
introducing a set of P and S waves for each P and each S wave introduced in 
the second stage and they are called the PPP and PPS, and PSP and PSS, ... 
for a total of eight. In particular the potentials of the SPP and SPS waves, 
® CDD and <t> , respectively, are introduced at z = 0 to cancel the stresses o' 

Dr r DrD 22 

and cr generated from the potential , corresponding to the reflection of 
rz SP 

the SP wave at a stress- free plane boundary at z=0 . 

The superpositions of the stresses generated from the potentials 
introduced in the first, second and third stages satisfy the boundary 
conditions at z = G only. The procedure of adding new reflected waves at each 
successive stage can be continued as above. At each stage, a new set of P and 
S waves for each P and each S wave introduced at the previous stage is added. 
This process is represented in Fig. 2 as a binary tree. At each stage, the 
boundary conditions at one face are satisfied and the boundary conditions at 
the other face are destroyed. However, as shown below, the waves introduced 
for increasingly large numbers of stages, namely the waves with very large 
numbers of reflections, vanish. Thus, the superpositions of all the stresses 
generated from the wave potentials from the first stage to infinity satisfy 


the boundary conditions at both z=0 and z=h simultaneously. Also, again as 
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shown below, the waves with increasingly higher numbers of reflections are 
observed at a receiving point with increasingly larger time delays? and thus 
considering a finite number of waves is sufficient for the exact analysis of a 
transient response up to a fixed time. 

In studying the wave propogation in a multi- layered elastic medium, 
Spencer Cl 41 analyzed multiple reflections and transmissions of axially 
symmetric waves at plane boundaries. Following the analysis given in C141 , in 
general the potential of a particular wave which has a unique chain of P's and 
S's can be shown to be Cl 31 


0 ) 


n / j 



E <€ ) Q(€ ) J ( 2tt€ ) 

x i 


h 

exp C 2tt i C < ( p-m ) 

x i 



+ <s-m „) J k 2 -S 2 ) ) +m — i b 2 ~S 2 3 } d€ (11) 

2 

where d>^ denotes the particular wave potential and is called the j-th wave 

v 

potential at the n-th stage. The unique chain for the particular wave 
potential is found by converting the decimal integer j to the corresponding n 
bit binary number and then assigning P's for 0's and S's for 1's in the binary 
number, n can take values from 1 to oo , and j can take values from 0 to 2 n -l. 

For example, the chain for the wave whose potential is represented by can 

be found as follows. First, j = 5 is converted to the n=4 bit binary number 
0101, and then P's for 0's and S's for l‘s are assigned in the binary number, 
which results in the PSPS wave. 

ECS ) in eqn. (11) is called the excitation function and is equal to 
Ep(C) or if is constructed from the half-space solution <t>p or 
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(that is, if the chain for started with a P or an S wave), respectively. 

Q(€ ) is called the product of the reflection coefficients and is given below. 

The numbers p and s represent the totajl number of P's and S's in the chain, 

respectively, and n=p+s. m^ is equal to 1 if the last reflection occurs at z=0 

(that is, if n is odd) and if is a P wave potential (that is if the chain 

ends with P) . Otherwise m is zero, m is equal to 1 if n is odd and if <t> 

1 2 n/j 

is an S wave potential (that is, if the chain ends with S) . Otherwise m is 

2 

zero, m is equal to +1 or -1 if the last reflection occurs at z=0 or z=h (that 

is, if n is odd or even), respectively, b is equal to 1 or k if 0 is a P or 

n/j 

S wave potential (that is, if the chain ends with P or S, or equivalently if j 
is even or odd), respectively. 

The product of reflection coefficients is given by 


- l 


* CQ pp <€)l 


Q oe (€ > ^ 
PS 


C iSQ e (€ ) 3 
SP 


C Q SS (* )3 (-1) (12) 

where j^, j , j , and j are the total number of P to P, P to S, S to P, and S 
to S reflections, respectively; and j and j are the total number of P to S 
and S to P reflections only at z=h, respectively. Q pp/ Q pg , Q S p' and Q ss are 
the P to P, P to S, S to P, and S to S reflection coefficients, respectively, 
and are given by 
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4€ 


2 J 1-* 2 J k 2 ~S 2 - C2€ 2 -fc 2 ) 2 


q pp (€) - Q ss (f > = 


G<€ ) 


( 13 ) 


4€ 


2 2 2 
sJl-5 <fc -2S ) 


°ps <5 ’ = 


G(€ ) 


( 14 ) 


- 4€ 


vfc 2 -S 2 ( 1c 2 -2S 2 > 


Q sp <s> = 


G(f ) 


( 15 ) 


j , j-/ j-/ j.# j c , and j can be found from the specific wave chain. For 
1 2 3 4 5 o 

enample for the PPSSPS wave, they are 1, 2, 1, 1, 1, and 0, respectively. Also 
for this particular wave E<£)=E _<§>; p=3 and s=3; m =0 and m =0; m=-l and b=k . 

So, as discussed above the superpositions of the stresses generated from 
the potentials introduced in the first stage as the half- space solution and 
all the possible reflected wave potentials introduced in the successive stages 
satisfy the boundary conditions at both z=0 and z=h as given by eqns. (1) and 


(2) , respectively . 
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3 FREQUENCY RESPONSE OF THE PLATE 


In general, a receiving transducer located on the plate observes a 

physical quantity such as stress, strain, or displacement generated from 

incident wave potentials. In this study, the SWF configuration is considered 

and it is assumed that a circular receiving transducer , which is separate from 

the transmitting transducer, is located at the top face of the plate. As in 

C73 , it is assumed that the receiving transducer produces an electrical 

voltage proportional to the normal stress, o' , generated from an incident 

zz 

wave potential averaged over its contact area. The receiving transducer 
therefore observes only the reflected waves at z=h which are introduced at the 
even numbered stages. 

It should be noted that the existence of the receiving transducer is not 
taken into consideration in Sec. 2, where stress waves in the plate are 
studied. Then, the question arises as to whether or not the receiving 
transducer affects the waves in the plate at the point of reception. The waves 
in the plate are affected by the presence of the receiving transducer after 
they reach the receiving transducer. Thus, incident waves on the receiving 
transducer are not affected, but waves that reflect at the receiving 
transducer contact area and continue to travel are affected because stress- 
free boundary conditions are not satisfied on the receiving transducer contact 
area. In this study, only the incident waves on the receiving transducer are 
considered and thus the results in Sec. 2 are used. 

Since a steady-state harmonic excitation having the magnitude of unity 
is considered and the plate is assumed to be linear and time invariant, the 
quantity observed by the receiving transducer may be called the frequency 
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response of the plate, or simply the frequency response. The quantity observed 
at a point on the receiving transducer contact area is called the frequency 
response at a point and is given by 

n =oo , j a J 

h . y H . 

n/ 3 

n = 2 , j = 0 

where H is the frequency response at a point due to the j-th wave at the n~ 
th stage. After dropping the steady- state term eHp(-iwt), H is equal to 

the normal stress cr^ generated from the wave potential & n/ j by the relations 
given in Table 1 and evaluated at the receiving point at z=0. It is given by 
C 1 31 


n=2 ,4,6 


J=2 n -1 


( 16 ) 


H , . = 2tt 
n / j 


-f 

X •'0 


ECS) QCS) VCS) J ( 2tt 5 ) exp 


r 

X, 


h 2 

C 2*rt i (ps| 1-5 


+ s \ k 2 -S 2 ) 


dS 


( 17 ) 


where VCS ) is called the response function and is equal to the response 

functions for the P wave V (5), or for the S wave V (5 ) if the potential O 

P S n h 

is a P or S wave potential, respectively, and where 


V p <€) = 25 2 - fc 2 


( 18 ) 



v *> 


(19) 



Only even values of n must be considered in eqn. <17) and all subsequent 
equations in this paper. 

The frequency response of the plate, H, or the frequency response of the 
plate due to the j— th wave at the n— th stage, can be found by averaging 

H or respectively, over the contact area of the receiving transducer. 

For a receiving transducer of radius a^ it can be shown that C131 


2 * t +a 

H =s I H . , r (x(r) dr (20) 

n / j ^ 2 J. n/j 

tt a _ ^ l - a 

2 2 


where 


-1 

<x ( r ) = tan 


C r 2 ( 21 2 +Za 2 2 -r 2 ) - ( « 2 
__ — ~ ~ 

r “* 2 + * 




( 21 ) 


and Z is the distance between the transmitting and receiving transducer axes. 

The integrals in eqns . (17) and (20) cannot be evaluated analytically in 
closed form. Thus they are evaluated approximately as discussed in the 


following sections . 
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4 ASYMPTOTIC EVALUATION 


The integral in eqn. (17) can be evaluated asymptotically by the method 
of stationary phase C153 . Accordingly, the following can be derived C133 


S E<! 0> Q({ 0 > V( V l9<4 o’ 

. ~ - 2 IT i 

i 

X 1 J (r/X 1 ) S Q I g " <£ Q > l 


<22 ) 


where € is the stationary point and which is defined explicitly below. g<€) 
is the phase function and is given by 


g ( € ) = 2tt <p >il-€ 2 + s J Ic 2 -€ 2 ) + 2-tt 5 (23) 

N x i 

I I denotes the “magnitude of“ and g“<€ ) is the second derivative of g<€ ) and 
is given by 


d 2 g -h 



2 3/2 

< 1 > 


2 2 3 /2 


< 24 ) 


The first derivative of the phase function vanishes at the stationary 
point, requiring that the stationary point satisfies the following nonlinear 
equation 



'I'-V » 


< 25 ) 
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For the interpretation of the equations above, the multiply reflected 
ray corresponding to the PSSP wave is constructed geometrically between the 
origin O and a receiving point M in Fig. 3 (a) as an example. A P ray 
emanating from the origin at an angle © 1 reflects as an S ray at , then as 

an S ray at N 2 and finally as a P ray at N g and then propagates to M. 
Reflection angles are determined by Snell's law C123 and thus all the P and S 
rays have reflection angles of ©^ and © 2 , respectively, and are related by 

sin ©^ = k sin ©2 < 26 > 

Using the ray geometry and eqn. ( 26 ), it can be shown that 

€ 0 = sin e i = k sin e 2 <27> 

For unmixed waves, namely for either p=0 or q = 0, ^ Q in eqn. (25) can be 
evaluated by substituting eqn. (27) into eqn. (25) and using trigonometric 
functions. Otherwise, numerical techniques such as the method of bisection 
C16] can be used. 

Further , it can be shown that the total distances travelled by the waves 
as P and S rays, R^ and R^, respectively, are given by 
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Ph 

R = <28) 



and the corresponding time delay is given by 


t _ + 

P' s „ « 


g(s 0 > 


( 30 ) 


Notice that the waves with common p and s, simply the waves with p,s, 

have the same time delay and the time delay appears as the linearity constant 

in the phase function which linearly depends on frequency. For example the 

waves with p=3, s=l are the PPFS, PPSF f FSFP, and SPPP waves. The frequency 

response at a point due to the waves with p,s, which is denoted by ^ , can 

be found simply by the superposition of the corresponding frequency responses. 

The frequency response of the plate due to the waves with p,s, H , is then 

P / s 

the average H over the contact area of the receiving transducer and can be 

P/S 

found similarly by eqn. (20). 

Also, it can be shown that the reflection coefficients given by eqns . 

(13), (14), and (15) are equal to the corresponding reflection coefficients 
for plane waves C123 when they are evaluated at the stationary point given by 
eqn. (27). This suggests that a receiving point observes an incident wave 


asymptotically as a plane wave propogating and reflecting at the bottom and 


top faces of the plate in the ray directions as constructed above. 
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The denominator function in eqn. (22) indicates geometrical attenuation 

for a given combination of r,h,p, and s. It can be visualized by considering a 

hypothetical multi- layered half-space as shown in Fig. 3 (b) . It is assumed 

that the layers which are identical to the original plate under study and are 

bonded together in a hypothetical way such that an incident wave in one layer 

produces no reflected waves but produces transmitted longitudinal and shear 

waves in the next layer and the transmission coefficients are the same as the 

reflection coefficients as if the incident wave reflects at a stress- free 

plane boundary. At the n-th layer in such a medium, only the waves introduced 

at the n-th stage in the binary tree in Fig. 2 could exsist . The top layer can 

be thought of as the original plate and each subsequent layer can be thought 

of as the image of the layer immediately above it . A receiving point M in the 

original plate and the images of it M^, M^, ... are shown in Fig. 3 (b) 

The PSSF ray between the origin and is also shown. So, the geometrical 

attenuation in eqn. (22) can be observed as due to the spreading of the waves 

with the distance from the source, and wave mode conversions during 

transmissions at the boundaries of layers. It can be shown that for unmixed 

waves, eqn. (22) represents a spherical wave whose magnitude is a function of 

directivity and decreases with the reflection coefficient Q (§ ) = Q (£ ) each 

FP 0 SS 0 

time it passes a layer, and the results obtained in C7] can be recovered. 

As observed from Fig. 3 (b) qualitatively, although some image points 
may be in the near field of the source, at sufficiently large distances from 
the source there are always image points in the far field. The value of the 
frequency response at a point as given in eqn. (22) is approximately valid for 
points in the far field. 

It should be noted here that in the limit as r approaches zero. 
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J ^(2*n€ a^ /X^ )/ nJ r€ q approaches *na ^/(X^\l nh) , and thus eqn . (22) is finite for 

r=0. Also, it can be shown that in the limit as n goes to infinity, H , . in 

n/j 

eqn. (22) goes to zero. 

4 . 1 Near Field and Far Field 

The analogy to the acoustics problem considered in C71 is extended to 
obtain a condition for the asymptotic evaluation in eqn. (22) to be a valid 
approximation to eqn. (17). The pertinent acoustics problem is that of a rigid 
piston vibrating in a rigid baffle. The exact formulation in integral form and 
the asymptotic solution for this problem are given in C173 . They are similar 
to eqn. (17) and (22), respectively, for unmixed waves. 

The exact solution for the radiation pressure along the piston axis can 
be obtained and is given in Cl 81 . The amplitude of the pressure along the axis 
(z axis in Fig. 1 (a)) fluctuates near the piston and eventually monotonically 
approaches the asymptotic solution. The near field and far field distances 
defined for such a problem are reviewed in C71 . The condition for the 
asymptotic solution to be valid is stated as 

zX„/a 2 > K ( a /X _ , E ) (31) 

0 F 0 a 

where z is the distance between the origin and the receiving point on the 
piston axis; is the wavelength in the acoustic medium; a is the radius of 

the piston; and is the far field distance constant which depends on the 
dimensionless parameters, a/X Q , and an acceptable error, E^ . is defined as 
the error between the asymptotic pressure amplitude P and the exact pressure 
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amplitude P relative to P 

ex ex 


E 

a 


( P - P ) 
as ex 


P 

e x 


< 32 ) 


Calculated values of for a range of its independent variables are 
given in Table 2. Notice that for fixed error, becomes constant for a/X Q 
much larger than l . Also the error decreases with increasing K for fixed 

a/ V 

The condition for eqn . (22) to be valid for unmixed waves can be stated 

2 

by analogy to the acoustics problem. zX Q /a in eqn. <31> is replaced by 
2 2 

° r R 2 X 2 /a i and the value of f ° r s P ecific values of or *1^2 

and E is taken from Table 2. X- is the S wavelength and is equal to X 4 /k. 
a 4 1 

For mixed waves, the analogy may be extended in a way, which also 
covers the condition for unmixed waves as special cases. Let the condition be 


‘Vl* R 2 X 2 ''*1 


K ( a /X ,E > 
F 1 av a 


< 33 ) 


where X ^ is the average wavelength defined by 


X 


a v 


( pX ^ + s X ^ ) / (p + s) 


< 34 ) 



21 


5 NUMERICAL EVALUATION 


The integral in eqn. (17) can also be evaluated by numerical methods 
[163. The integral is in the complex € -plane because the integrand has 
singularities at the branch points 5=1 and 5 =k, and at the Rayleigh pole 5=5 

R 

on the integration path from 0 to w C193 . The pole is of order n=p+s because 
all the reflection functions given by eqns. (13), (14), and (15) have the 

Rayleigh pole of order 1 as well as the excitation functions given by eqns. 

(6) and (7) . The introduction of appropriate branch cuts and indentations of 
the integration path on the real axis near the singular points is discussed in 
C83 . Consistent with the branch cuts introduced, the principal values of J 1-5 ^ 

I 2 2 

and \J k -5 should be considered. Thus the exponential function in the 
integrand in eqn. (17) is purely oscillatory for 0<5<1; the part associated 

I 2 I 2 2 

with v) 1-5 monotonically decays, and the part associated with nJ k -5 is 
oscillatory for 1<5 <k; and monotonically decays for 5>k. 

The integration interval from 5=0 to oo is divided into subintervals and 
as indicated in Table 3, either the trapezoidal rule or Gauss quadrature C163 
is used for each subinterval. The integration is stopped at some large value 
oi 5 =5 k , where 5^>k/ which is justified by the fact that the integrand is 
exponentially decaying for 5>k. The trapezoidal rule is chosen because it is 
the simplest method which works well for oscillatory integrands Cl 63 . The 
Gauss quadrature is chosen in the vicinity of the singular points because it 
does not require the evaluation of the integrand at the end points of the 
interval, thus avoiding evaluations at the singular points. 

The convergence and stability of the procedures above depend on the 


exponential decay for 5>k. As shown below, if the decay is not sufficiently 
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rapid, the effect of the Rayleigh pole is observed as an instability. 

Increasing h/X^ , p, or s increases the rate of the exponential decay, thus 

improving the convergence and stability; however, increasing p or s also 

increases the power of the pole, thus promoting instability. Thus, there is a 

compromise associated with increasing p or s in order to achieve a dominance 

of the exponential decay over the algebraic instability. For a given set of p 

and s, there is a minimum h/X , (h/X ) , below which the procedure above 

l l mm 

gives unstable results. By examining the exponential decay and using the fact 

J 2 I 2 2 

€ -1 > \1€ -1c for € >fc, it can be shown that for the cases with fixed 

n=p+s the case with p = 0, s=n gives the largest <h/X ) 

1 min 

The integration over the receiving transducer contact area as given in 
eqn . (20) can also be performed numerically using the trapezoidal rule . 
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6 EFFECT OF LOSSES 

So far, it has been assumed that the plate material is perfectly elastic 
and all energy is conserved during wave propogation as only geometrical 
attenuation due to spreading and wave mode conversions is present. However, 
wave propogation in real materials is never entirely conservative as a certain 
amount of energy loss always occurs during propagation and reflections. The 
mechanisms of these losses are discussed in C203 . In this study, it is assumed 
that attenuation is small and the effect of the attenuation can be introduced 
by multiplying eqn. (17) by expC-Ccx^* o^)] where and o< 2 are the 

frequency dependent attenuation constants for P and S waves, respectively 
C203 . oc^ and c* 2 are material properties and have the units of nepers per 
centimeter C63 . 
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7 RESULTS AND DISCUSSIONS 

The frequency response of the plate in the SWF configuration is 
formulated as the superposition of the frequency responses due to reflected 
fields. The response of the plate to an arbitrary pulse is thus the the 
superposition of the responses due to reflected fields to that same pulse, 
simply reflections, and can be found by Fourier analysis C21I . 

The frequency responses due to the fields with common p,s have the same 
phase which depends linearly on frequency; and thus the corresponding 
reflections whose superposition is called the reflection with p,s have the 
same time delay in the transient response. If the plate is sufficiently thick, 
then the time difference between the time delays of the significant 
reflections with different p,s may be so large that each significant 
reflection can be observed in the output separately as in C71 . When the plate 
becomes thinner, the reflections overlap and it is no longer possible to 
obtain information about an individual reflection directly in the output. 
However, it may be possible to do so by using homomorphic signal analysis 
[13,22:1. 

Information about one reflection is useful because it is a prerequisite 
to analyzing the entire output quantitatively and also because it may be far 
easier to interpret the changes in individual reflections due to defect states 
of the plate. This is especially true if the reflection is in the far field, 
where a receiving point observes the incident field as a plane wave. Although 
the thickness of the plate may be small, reflections with large enough p or s 
satisfy the far field condition given by eqn . (34). So, information in the far 
field may be obtained from some portion of the output signal. 
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The frequency response at a point due to the field with p f s H depends 

P/S 

on the dimensionless parameters p, s, k, cx \ cx X , a /X , h/X , and r/h . The 

11 is is 11 1 

frequency response of the plate due to the field with p,s H also depends 

P • s 

on the dimensionless parameters p, s, k, cx^Xj , 1 *1^1' and 

a /a . Unless otherwise stated, it is assumed that the plate material is 

^ i 

lossless, namely cx^ =0x^=0, in the following calculations. Also, aluminum for 

which k=2.02 is considered as the plate material. 

The behavior of IH I , IH I , and IH I versus h/X is plotted in 

* / w 6,0 3,1 1 

Figs. A <a), (b) , and (c) , respectively, for the values of a^/X^ of 3 for all 
cases and for r/h of 0 , 0, and 2.5, respectively. The magnitudes calculated 
using the asymptotic formula in eqn. (23) are indicated by the filled squares. 

The magnitudes calculated by evaluating the integral in eqn. (18) by the 
numerical procedures described in Sec. 5 are indicated by solid curves. 


The asymptotic magnitude in Fig. A (a) shows spherical wave behavior. 


and thus it starts to “blow up" for small h/X^ , namely, near to the source. 

The numerical magnitude has approximately the same value as the asymptotic 

magnitude for large h/X^ (for h/X^>14.4 with the criterion given by eqn. (3 A)) 

and it is oscillatory for smaller h/X The numerical magnitude also "blows 

up" at (h/X ) . ~0.48 because of the effect of the Rayleigh pole as discussed 
1 min 

in Sec. 5. Similar discussions can be given for Figs. A (b) and (c) . 


As discussed in Sec. 5 (h/X ) . is controlled by p and s. For the cases 

1 min 

with fixed n=p+s, the case with p=0, s=n gives the largest (h/X ) . . 

1 min 

(h/X ) . has been evaluated for the cases with increasing s and fixed p=0, 
l mm 

a^/X^S, r/(sh) = 0.7. It increases with increasing s and gives a peak value at 
aproximately 1.6 for s-10 and then decreases as suggested by the discussions 
on the stability of the numerical procedure in Sec. 5. The critical h/X^ , 
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<h/\j) cr , which is the largest (h/X^)^.^ among all the cases, is approximately 
1.6. (h/X^)^ limits the applicability of the formulations in this study. For 

h/X^ <(h/X^)^ , expanding the plate modes derived by Rayleigh- Lamb frequency 
spectrum, which is more suitable for low frequencies, may be used. 

Some parametrized plots which determine the particular field whose 
frequency response has the maximum magnitude compared with other reflected 
fields are given in Figs. 5-8. These plots may also be used in the spectral 
analysis of transient signals to estimate the particular reflection which has 
the maximum magnitude compared with other reflections at a particular 
frequency. The reflection with the maximum magnitude is important because it 
has the maximum signal to noise ratio in an experiment . 

Since the transmitting and receiving transducers considered in this 
study are longitudinal, the reflections with only F waves, namely with s=0, 
will be dominant in the output [7,133; and thus only the fields with s = 0 are 
considered. Again, k=2.02 is used, but the results can be generated for other 
materials . 

In Fig. 5 <a) , a point receiver, namely a /a =0; a lossless material, 

A 1 

namely o^X^O; and a 1 7X^=1 are considered. The abscissa and the ordinate are 
h/X^ and £/a^, respectively. Then all the parameters in this analysis except p 
are fixed for each point in the <h/X 1# */a ) plane. 

Suppose the magnitude of the frequency response at a particular point 

due to the field with p = p is maximum among all the possible fields with p = 2 

max c 

t0 ** * P max be found *° be the saTne for each point in a region bounded by 

solid curves in Fig. 5 (a) (and is indicated in each region. For example, 

=1 for h/X. =6, */a =5. 
max 1 1 

In generating Figs. 5 8, discrete points where h/X^ and ^/a^ change from 
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2 to 40 and 2 to 10, respectively, with increments of 1 have been considered. 

Then IH I (actually 1H _l for a point receiver) has been calculated for 

each point by changing p from 2 to a sufficiently large number, in increments 

of 2. The asymptotic formula given by eqn. (23) has been used if the condition 

given by eqn. (34) has been satisfied with the value of K for E =5% in Table 

F a 

2. The numerical procedure discussed in Sec. 5 has been used otherwise. 

IH I may have diverged with increasing p at the beginning, but it has 

p ( o 

finally converged to zero due to geometrical attenuation as discussed in Sec. 

4. Thus, the sufficiently large number for p has been set by checking the 
convergence and when 1 Q I has become negligible (20 dB down) compared with 
the maximum magnitude. The curves dividing the regions having the same p 

max 

have been then generated by linear interpolation . 

Similar plots are given in Figs. 5 (b) and (c) where only a^/X^ as 3 and 

5, respectively, is changed. For the plots in Fig. 6, again lossless material 

but a^/a^ =1 are considered, and a^/X^ is 1,3, and 5 for Figs. 6 (a), (b) , and 
(c), respectively. In Fig. 7 the conditions are the same as in Fig. 5 <b> 
except <*^X^ is changed to 0.005 and 0.01 in Figs. 7 (a), and (b) , 
respectively. Similarly, the conditions in Fig. 8 are the same as in Fig. 6 
(b) except is changed to 0.005 and 0.01 in Figs. 8 (a) and <b> , 

respectively . 

There are several observations which can be made from these plots. 

1. In the plots in Fig. 5, ail the fields with p=p are in the far 

max 

field, p decreases with increasing h/X 4 and increases with increasing */a 4 
max 1 1 

and a^/X^ . As a^/X^ becomes larger or h/X^ becomes smaller, changes in h/X^ or 

produce larger changes in p 
i max 

In order to understand these and the following observations, consider 
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the hypothetical multi- layered half space in Fig. 3 <b) again, and consider 
only the P waves in the medium. Then the field with only P waves observed at 
an image point M is the field with p,s = 0 observed at the receiving point M in 
the original plate as discussed in Sec. 4. Since the excitation with a^/X 1 is 
longitudinal in the z direction, the field with P waves only dies out at the 
points located away from the source. As discussed in Sec. 4, in the far field 
it is a spherical wave whose amplitude is a function of directivity and 
decreases due to mode conversions each time it passes a layer. The polar 
diagrams of the directivity function for several a^/X^ are given in C73 . It is 
mainly controlled by the function f(u>= I J^uJ/u I where u=2*n€ Q /X ^ and 5 Q =sin 
©^ . The properties of f(u> is discussed in E233 . For sufficiently large a^/X^, 
the directivity function has a main lobe where it has a maximum at ©^=0 and 
decreases to zero, and beyond which it has successive side lobes where it 


reaches a maximum and decreases to zero with increasing ©^ . The maximums in 
the side lobes decrease with increasing angle ©^ and they are much smaller 
than the maximum at ©^=0. The number of side lobes increases with increasing 
a^/X^ and the main lobe is confined to a smaller angle. Thus, the decrease in 
the directivity function with increasing angle becomes faster. For small 


a^/Xi , the directivity function has only a main lobe and it never becomes 


Now notice that an image point M in Fig. 3 <b) is located at a larger 

P 

distance R but at a smaller angle ©, with increasing p. Then I H _1 tends to 
1 1 P / u 

decrease due to larger geometrical attenuation because of the larger distance 
and larger number of mode conversions, but also it tends to increase due to 
the directivity function because of the smaller angle with increasing p. So, 


there is a compromise in the value of p 

max 
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As h/Xj^ increases or Jf/a^ decreases © 1 becomes smaller for smaller p and 

P tends to decrease because of the directivity function. As a./X, increases 
max 1 1 

the value of the directivity function for larger ©^ , namely for smaller p, 

becomes smaller and p tends to increase. For smaller h/X 4 , changes in h/X„ 

max 1 1 

or €/a^ produce larger changes in ©^ and thus in the value of the directivity 

function and p Also, for larger a. /X,, changes in h/X« or If a. which 

max 11 1 1 

accompany changes in ©^ produce larger changes in the directivity function and 
thus in p 

max 


Let 

the fields with p,s=0 

be in the 

far 

field 

for p>p^ From 

far 

the 

condition 

given 

in eqn. (34) 

it can 

be 

shown 

that p^ increases 
far 

with 

increasing 

vs 

and decreasing 

h/X ^ or 

€/ v 

Also 

as discussed above 

P 

max 

increases 

with 

increasing a^/X^ and 

< /a i 

, but 

with decreasing h/X 

t . The 


results show that always p maK >P£ ar for range of parameters considered in 

Fig. 5. Notice that for decreasing */a . p decreases and p^ increases. 

1 max far 

Thus it is expected that for sufficiently small €/a . p <p^ . But this does 

1 max far 

not happen for ^/a^2. 

It can also be shown that © for the fields with p , which is denoted 

1 max 

by (© ) , is in the main lobe of the directivity function for all the cases 

1 max 3 

considered in Fig. 5. 

2. Similar observations as for Fig. 5 can be made for Fig. 6. The 

conditions for Figs. 6 (a) , (b) , and <c> are the same as for Fig. 5 (a), (b) , 

and (c) , respectively, except the receiver is such that a /a =1 instead of a 

point receiver. It is observed from Figs. 5 and 6 that p tends to increase 

max 

for a^/a ^ = 1 compared with a^a^O. 

For a^/a^X), the complex frequency response ^ is integrated over the 
receiving transducer contact area to find H „ and then I H ^ I . Suppose 

PiO p, 0 
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IH „ I is the average of IH -I over the receiving transducer 
P,0 p,0 

contact area. Also, suppose that IH „ l is the value of I H _l evaluated at 

p,0 0 p,0 


the center of the receiving transducer. For IH _l ~IH _l n , which can 

P , 0 P ; U U 

be justified for the fields with © in the main lobe, IH «l 5 IH _! n if 

1 p, 0 p , 0 0 

there is a change in the phase of H over the averaging area . It can be 

p # 0 

shown that for smaller p the phase changes more over the receiving transducer 
contact area by considering that the imaginary point in M in Fig . 3 <b) is 

p 

located at a smaller distance and at a larger angle from the source, and thus 

IH I decreases more compared with IH n l . So, p tends to increase for 
p , 0 p , u max 

V*1= 1 - 

3. The conditions for Fig. 7 are the same as for Fig. 5 (b) except <*^X^ 

is 0.005 and 0.01 in Figs. 7 <a> and (b) , respectively , instead of lossless 

material. As observed from these plots, p decreases with increasing X . 

max 11 

The plots with attenuation show irregularities for large h/X^ and {/a^, where 
p c * oes no * necessarily decrease with increasing h/X^ or decreasing £/a^ as 


in lossless material. The regions with irregularities move towards smaller 

h /X, or */a 4 with increasing o< X 4 . 

1 1 11 

As assumed in Sec. 6, I H^ decreases exponentially with the product of 

the attenuation constant and the distance of the imaginary point M from the 

P 

origin in Fig. 3 <b) . Thus, IH I with larger p decreases more with 

p # 0 

increasing cx x . So, p tends to become smaller with increasing X. . 

1 1 max 1 1 

Although (0 . ) for the field with p is in the main lobe of the 
1 max max 

directivity function for in Fig. 5 (b> , it may jump to the side lobes 

for sufficiently large h/X^ , €/X^, and 0< 1 X 1 , and irregularities appear as in 
Fig. 7. While attenuation pushes p down, namely (0.) up, (0 4 ) 

approaches to the critical angles where the directivity function is exactly 
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zero. Near these angles there are sharp changes in l H I and thus p jumps 

p,0 man 

to bigger or smaller values . 

4. The conditions for Fig. 6 <b) f 8 (a) and (b) are the same as for Fig. 

5 (b) , 7 (a) and <b) , respectively, except a^/a^sl instead of a point 
receiver . 

For a 2 /a i = l' p max decreases with increasing again but the 

irregularities appeared for a point receiver do not exsist any more. This is 
because I R | is not exactly zero near the critical angles due to the 

p,0 

averaging over the receiving transducer contact area and there are no sharp 
changes . 

5. Overall, Figs. 5-8 suggest that larger h/X^ , smaller a^fX^, and for 
attenuating material larger a /a increase the stability of the output in the 
SWF configuration in terms of the reflection with the maximum magnitude. The 
stability means small changes in parameters produce small changes in the 


output . 



32 


8 CONCLUSIONS AND RECOMMENDATIONS 


Steady-state harmonic stress waves in an isotropic elastic plate excited 
on one face by a circular transducer were analyzed theoretically. It was 
assumed that the transmitting transducer transforms an electrical voltage into 
a uniform normal stress at the top face of the plate. First, the radiation 
into a half- space, which introduced longitudinal (P) and shear (S) waves into 
the plate, and then their successive reflections at the bottom and top faces 
of the plate were considered to satisfy the boundary conditions. 

A separate circular receiving transducer, also located at the top face 
of the plate, was considered. It was assumed that the receiving transducer 
produced an electrical voltage proportional to the average spatially 

integrated normal stress over its face due to an incident wave. 

The frequency response at a receiving point due to a multiply reflected 
wave was formulated in integral form and its asymptotic behavior was given. 
The far field condition for the asymptotic solution to be valid was also 
discussed. This condition suggested that although the thickness of the plate 
may be small, the waves which reflect sufficiently many times satisfy the far 
field condition. 

A numerical procedure was given to evaluate the frequency response at a 
receiving point due to a multiply reflected wave in the near field. Its 
stability and convergence were discussed. Also, exponential decay was 

introduced to account the material attenuation. 

Calculations were done for aluminum plates. It was found that the 
numerical procedure becomes unstable for h/X^<1.6 where h is the thickness of 
the plate and X^ is the P wavelength. 
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Parametrized plots which determine the particular wave whose frequency 
response has maximum magnitude compared with other multiply reflected waves 
were given for a range of values of dimensionless parameters in the analysis. 
The effects of changes in the values of the parameters were discussed. 

This study is part of an overall effort to develop quantitative analyses 
of ultrasonic nondestructive evaluation parameters such as the stress wave 
factor (SVF) . The spectral analysis of the SVF signal should benefit from this 
study. The parametrized plots given in this study or generated using the 
results of this study may be used to estimate some parameters such as the 
material attenuation or the effective radius of the transmitting transducer 
from the spectral analysis of the SVF signal in terms of individual 
reflections . 

Finally, stress wave transmission characteristics of plates in the SVF 
configuration should be analyzed using the modal analysis (Rayleigh- Lamb 
frequency spectrum) for small values of h/X^ where the numerical procedure 
given in this study becomes unstable. 
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TABLE 1 Displacements* and Stresses* in Isotropic Elastic Solid with Axial 
Symmetry and Zero Rotational Displacement C123 . 
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TABLE 2 


Calculated Values of for Analog Acoustics 
Problem. 
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TABLE 3 Subintervals and Corresponding Numerical Integration 
Methods Used to Evaluate the Integral in Eqn. (18). 


Subinterval 

Numerical 



Integration 

From 

To 

Method 

0 

l-«* 

b 

Trapezoidal 

l-€ 

l 

Gauss C 

1 

u« 

H 

l+« 

k-« 

Trapezoidal 

lc-i 

k 

Gauss 

k 

< 

it 


R 


< 


<• 

R 

1 R d 


<i 

2 1 me 

Trapezoidal 

i. 


M 

2 

3 2 me 

If 


; 

II 

II 

r 4 

f ml 4 

II 

k-t 

k k-l inc 


S* 

0* 

Naglectad 


a «<<l f typically «*0.01. 

b Number of steps is increased by powers of 2 and the 
convergence is checked. When the integral converges 
to a value with an acceptable error, typically to 
within of 0.001, the integration is stopped. 

c Typically 8 point Gauss-Legendre quadrature is used. 

d Typically, i. =1. 

me 

• The integration is stopped it 4 fc when the integral 
over the interval from S ^ to t becomes negligibly 
small, typically smaller than 0.001. 
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(a) (b) 


Fiq. 1 Schematics of (a) plate with axially symmetric excitation and 
(b) coordinate systems. 


Boundary Conditions 

Stage Waves Satisfied Only at 



A A A A 

3 PSPSPSPS z =0 

A A A A A A A A 


Fig. 2 Representation o< introducing reflected waves at successive 


stages as binary tree. 




Fig . 3 Schematics of 
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(a) 


sin d | -k Sin02 



(b) 


t> put. and th. PSSP ray and (b) hypothetical 
f-spac. and th. PSSP rav. 
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(C) X | 


Fig. 3 Parametrized plots for k=2.02. o^x^O, s»0. a^a^O (point receiver) 
and <*> a^/X^sl. (b) ^7X^3, and <c) a^/X^ssS. 
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sidered. It is assumed that the receiving transducer produces an electrical voltage proportional to 
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behavior of the frequency response at a receiving point due to a multiply reflected wave is given. It 
is found that a receiving point observes an incident wave asymptotical ly as a plane wave which prop- 
agates and reflects in the direction of the multiply reflected ray constructed geometrically using 
Snell's law. The far field condition, for which the asymptotic solution is valid, is also discussed. 
This condition suggests that although the thickness of the plate may be small, the waves which reflect 
sufficiently many times do satisfy the far field condition. A numerical procedure is given to evaluate 
the frequency response at a receiving point due to a multiply reflected wave in the near field. Its 
stability and convergence are discussed. Also, exponential decay is introduced to account for material 
attenuation. Calculations are done for aluminum plates. It is found that the numerical procedure 
becomes unstable for hA i <1.6 where h is the thickness of the plate and X -j is the P wavelength. 
Parameterized plots which determine the particular wave whose frequency response has maximum magnitude 
compared with other multiple reflected waves are given for a range of values of dimensionless parame- 
ters. The effects of changes in the values of the parameters are discussed. 
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